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ABSTRACT 

We combine spatially resolved ASCA temperature data with ROSAT imaging data to constrain the 
total mass distribution in the cluster A401, assuming that the cluster is in hydrostatic equilibrium, but 
without the assumption of gas isothermality. We obtain a total mass within the X-ray core (290 h^Q kpc) 
of 1.2^0 5 ^ 10^'' h^o ^^0 ^^ ^^^ 90 % confidence level, 1.3 times larger than the isothermal estimate. The 
total mass within r5oo (1.7 h'^Q Mpc) is M500 = O-Q^oJ x 10^^ ^^0^ M© at 90 % confidence, in agreement 
with the optical virial mass estimate, and 1.2 times smaller than the isothermal estimate. Our Afsoo 
value is 1.7 times smaller than that estimated using the mass-temperature scaling law predicted by 
simulations. The best fit dark matter density profile scales as r~^-^ at large radii, which is consistent 
with the Navarro, Frenk & White (NEW) "universal profile" as well as the King profile of the galaxy 
density in A401. From the imaging data, the gas density profile is shallower than the dark matter profile, 
scaling as r^^'^ at large radii, leading to a monotonically increasing gas mass fraction with radius. Within 



rsoo the gas mass fraction reaches a value of fgas = 0.21 



+0.06 



-0.05 "50 



h, 



-3/2 



(90 % confidence errors). Assuming 



that fgas (plus an estimate of the stellar mass) is the universal value of the baryon fraction, we estimate 
the 90 % confidence upper limit of the cosmological matter density to be ilm < 0.31, in conflict with an 
Einstein-deSitter universe. Even though the NEW dark matter density profile is statistically consistent 
with the temperature data, its central temperature cusp would lead to convective instability at the center, 
because the gas density does not have a corresponding peak. One way to reconcile a cusp-shaped total 
mass profile with the observed gas density profile, regardless of the temperature data, is to introduce a 
significant non-thermal pressure in the center. Such a pressure must satisfy the hydrostatic equilibrium 
condition without inducing turbulence. Alternately, significant mass drop-out from the cooling flow 
would make the temperature less peaked, and the NEW profile acceptable. However, the quality of data 
is not adequate to test this possibility. 

Subject headings: cosmology: observations - dark matter - galaxies: clusters: individual (A401) - 
intergalactic medium - X-rays: galaxies 



1. INTRODUCTION 

Determining mass components of clusters of galaxies is 
an important task in observational cosmology, since clus- 
ters form through the collapse of a large volume of pri- 
mordial matter, and as such may provide a representa- 
tive s ample of the universe as a whole (e.g. White et al 
1993). Assuming hydrostatic equilibrium, the total mass 
of a cluster can be determined from t he intracluster gas 
temp e rature and dens ity distributions ( Bahcall fc Sarazir] 



isothermal estimates (e.g. White & Eabian 1995). In this 
paper, we estimate the total mass for the A401 cluster, 
using the actual temperature pr ofile . Our method is es- 



1977, Mathews 1978). Until recently, most hydrostatic 



X-ray mass estimates have been made assuming that the 
gas is isothermal at the average broad beam temperature. 
ASCA observations provide spatially resolved X-ray spec- 
troscopic measurements of hot clusters, and yield the 2D 
temperature structure of clusters. Indeed, a large number 
of ASCA clusters shows that the t emperature declines wit h 
increasing radius from the center ( Markevitch et al. 19981 ), 
in qualitative accordance with hydrodynamic cluster sim - 



ulations (e.g. Evrard et al. 1996, Bryan & Norman 1997) 



This implies that the real hydrostatic mass at large clus- 
ter radii is smaller than that derived assuming isothermal- 
ity. Conooquently, the gao maoo fraction io larger and the 



'baryon catastrophe" even more pronounced , compared to 



sentially the one used for A2163 (Markevitch et al. 1996 ) 
and A2256 ( [Markevitch fc Vikhlinin 1997^ . We assume 
that the cluster is in hydrostatic equilibrium and model 
the dark matter component using several different func- 
tional forms. We use the ROSAT data to fix the gas den- 
sity profile and fit hydrostatic temperature models, as a 
function of dark matter density parameters, to the ASCA 
data. 

A401 (z = 0.0748) is suitable for measuring the dark and 
total mass distributions, since it is fairly bright and hot (~ 
8 keV) allowing accurate temperature determinations with 
ASCA. Also, the ASCA field of view covers the cluster to 
rsoo (the radius where the mean interior density equals 500 
times the critical density, approximately the radius inside 
which hydrostatic equilibrium holds). A401 also has been 
observed with the ROSAT PSPC, which gives an accurate 
estimate of the gas content of A401 and shows no obvious 
substructure and no significant deviations fro m azimuthal 
symmetry. The ASCA 2D temperature map (Markevitch 
et al. 1998|) shows no strong asymmetric variation. Thus 
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the assumption of hydrostatic equihbrium is likely to be 
valid. A401 is peculiar in the sense that it shows no signif- 
icant evidence of a cooling flow (see Peres et al. 1998, who 
find only an upper limit for the mass deposition rate) , even 
though it has a prominent cD galaxy at its center and no 
evidence of recent merger activity either in the tempera- 
ture map or in the X-ray image. However, the ROSAT HRI 
image shows a linear structure in the neighboring cluster 
A399, pointing t owards A401, as pos sible evidence of some 
past interaction (Fabian et al. 1997). The lack of a signif- 
icant cooling flow simplifies our analysis since we need not 
consider multicomponent temperature models in the cen- 
ter. We use Ho = 50ft.5okms^^ Mpc~^, fi = 1 and report 
90% confidence intervals throughout the paper. 

2. ROSAT ANALYSIS 

2.1. Data reduction 

The ROSAT data consist of two PSPC pointings of 
A401, rp800235n00 and rp800182n00, the former taken on 
July 29, 1992 and the latter on January 22, 1992. The 
total exposure times are 7457 s and 6735 s respectively. 
Reductions were carrie d out using Snowden's Soft X-Ray 
Background programs (Snowden ct al. 1994), which re- 
duced the total exposure to 11.7 ks. The spatial analysis 
was restricted to the band of 0.44 - 2.04 keV (Snowden's 
bands R4-R 7) to improve sensit ivity over the X-ray back- 
ground (see David et al. 1995). The non-cosmic X-ray 



background was subtracted using Snowden's code. All de- 
tector and telescope effects, including vignetting, the mir- 
ror support structure shadows and varying detector quan- 
tum efficiency, together with the aspect and livetime infor- 
mation are incorporated in exposure maps for each band. 
All analysis was done in each of the four bands, for both 
pointings, the resulting images were divided by the cor- 
responding exposure maps and combined. The surface 
brightness contour map (smoothed by a Gaussian with 
a = 1') is shown in Figure 1. The data show no obvi- 
ous substructure and no strong deviations from azimuthal 
symmetry, implying that the assumption of hydrostatic 
equilibrium is likely to be valid. 

2.2. Spatial analysis 

With our spatial analysis, we wish to address two ques- 
tions. First, how far from the cluster center can we sig- 
nificantly detect the cluster emission, and second, how is 
the cluster gas distributed. Figure |l| shows that the sur- 
face brightness (cluster -I- background) reaches a constant 
level at a distance of about 30' from the center {— 3.4 
h'^Q Mpc) so we measure the Cosmic X-Ray Background, 
plus any residual detector background, from the same im- 
age at radii between 30' and 52'. To exclude any con- 
tribution from the nearby cluster A399 we excluded an 
azimuthal sector centered at A401, directed at A399 (154° 
clockwise from the north), from our analysis. To determine 
the region to be excluded, we increased the angular extent 
of the excluded sector until the background level reached 
a minimum value, with a sector width of 160°. Thus, the 
contribution of A399 to the 30' - 52' annulus (and inner 
annuli as well) becomes negligible outside this sector. We 
use the minimum value estimated above, 1.8 x 10""* counts 
s~^ arcmin"^, as the total background value. 

We excluded point sources and background fluctuations 



from our analysis. Around a radius of 20' (18' - 27') it 
is not clear whether the small fluctuations are associated 
with cluster emission or not. The detector support rib near 
20' may be causing some residual effects. This led us to 
two schemes, where we masked out, in addition to unam- 
biguous point sources at other radii, i) all point-like sources 
between 18' - 27' (a conservative scheme) and ii) none of 
those sources (a non-conservative scheme), the truth being 
somewhere between. 

We generated a radial surface brightness profile in con- 
centric annuli ranging from 15" at the center to 10' at a 
radial distance of 50'. The signal-to-noise ratios as a func- 
tion of radius include 5% of the background value as a 
systematic error added in quadrature. In the conservative 
scheme, we can detect the cluster gas emission with 2.6 a 
significance at r = 25', and in the non-conservative scheme 
with 3.8 (T significance at the same radius, so we conclude 
that we detect the cluster gas at ~ 3a up to 25' (2.9 h^Q^ 
Mpc). In our further analysis we consider only the profile 
of the conservative scheme. 

We fitted the observed profile with the (3- model plus 
background 



7(6) = /( 




(-3/3+i) 



BGD 



(1) 



(Cavaliere & Fusco-Femiano 1976), where b is the pro- 
jected radius and background BGD is fixed to the value 
found above. We used XSPEC to convolve the surface 
brightness model through a spatial response matrix (con- 
structed from the ROSAT PSF at 1 keV, for a spherically 
symmetric on-axis source) and to compare the convolved 
profile with the data. The brightness in the two inner- 
most bins (r < 30") rises 2.1 a and 1.2 a, or 47% and 
12%, above the best fit /? - model, respectively. This be- 
haviour is consistent with Peres et al. (1998) who found 
that the A401 data, consistent with no cooling flow, did 
allow a mass accretion rate up to 120 Mq/ji in the center 
of A401 within a cooling radius of 0.7~!^q'j arcmin. Even 
though the central excess in our data is statistically not 
very significant, we excluded our two innermost bins from 
the /3 fit to prevent any bias towards small values of the 
core radius. We find an acceptable fit in the radial range 
0.5'-52' (see Figure |2] and Table |^), with best fit parame- 
ters and 90 % confidence errors of a^ — 2.56 ± 0.14 arcmin 
(= 294±16 h^fl kpc), P = 0.70±0.02, Iq = 5.6±0.3 x (lO'^ 
counts s~^ arcmin"^) at r = 0', with x^ = 49.7 for 55 de- 
grees of freedom. The confidence contours for Ux and /3 are 
shown in Figure y. Our values of a^, and /? are consistent 
with the results of another study of the ROSAT data of 
A401 (Vikhlinin et al. 1998). 

If we assume that the intracluster gas is isothermal and 
spherically symmetric, the best-fit parameters ax and /3 
will determine the shape of the gas density profile by the 
equation: 



Pgas{r) = Pgas{^) 




(2) 



A401 is hot and even the temperature variations such as 
those detected by ASCA do not significantly affect the 
brightness in the ROSAT band. 
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We obtained the normalization (as in Vikhlinin et al 
l998Upgas{0) = 1-6 X 10"Mq Mpc-3, or 1.1 xlQ-^" g 
cin~ , by equating the emission measure calculated from 
the above equation, with an observed value of 16.7 xlO^^ 
cm~^ inside a cylinder with r = 2 h^Q Mpc radius, cen- 
tered at the cluster brightness peak. 

The gas mass in this best fit model inside the sphere of 
rsoo (=15.1' = 1.7 h'i^Q Mpc, calculated in Section 5.1) is 



Mgas{< 7-500) = 2.01 ± 0.08 X 10^^ h 



14,-5/2 



50 



Ma 



(3) 



White & Fabian (1995) give a gas mass value of 1.32 ± 
0.07 X 10^"* Mq inside 1.3 h'^Q Mpc whereas our corre- 
sponding value at the same radius is 1.4 x 10^^ M©, con- 
sistent with theirs. The error interval in (0) corresponds 
to the 90% confidence region in a^ - (3 - space in our global 
fit (see Figure 0) . The gas mass error within rsoo is negli- 
gible with respect to the other errors in the quantities we 
are interested in this work. Therefore we ignore the above 
error in our further analysis. 

3. TEMPERATURE DATA 

3.1. ASCA 

The gas temperature distribution is obtained from the 
ASCA spectral data of A401, excluding a sector towards 
A399, as described in detail in Markevitch et al. (1998). 
The data are divided into four concentric radial bins, 0'-2'- 
5'-9'- 16' (0-0.23-0.57-1.0-1.8 h^p^ Mpc). The temperature 



approximately constant temperature profile. From the de- 
scription in their paper, it appears that they did not prop- 
erly include PSF scattering effects. This would have an ef- 
fect of diminishing the radial temperature decline, consis- 
tent with the difference of the two results. ROSAT PSPC 
data on A401 in the 0.2-2 keV band were also analyzed 
by Irwin et al. (1999), who derive a temperature increase 
with radius. Taking into account the PSPC calibration 
uncertainty which is very significant for determining tem- 
peratures of hot clusters with PSPC (see, e.g., Markevitch 
& Vikhlinin 1997a), their results would probably be con- 
sistent with ours. Among other clusters in the Markevitch 
et al. (1998) sample, ASCA data for A4059 and MKW3s 
were recently independently analyzed by Kikuchi et al. 
(1999) using an independent method. Although these au- 
thors do not find radial temperature gradients as strong 
as in Markevitch et al. (1998), they also do not detect 
the strong cooling flows that are known to exist in those 
clusters (e.g., Peres et al. 1998), which indicates possi- 
ble problems with their results. For another cluster from 
that sample, Hydra-A, the Markevitch et al. (1998) de- 
clining temperature profile is in good agreement with an 
analysis by Ikebe et al. (1996) using a different technique. 
A strong radial temperature decline was also recently de- 
rived for A2218 by Cannon, Ponman, & Hobbs (1999) us- 
ing yet another method. Still another technique developed 
by Churazov et al. (1996) yields a temperature structure 
in Coma and A1367 consistent with the Markevitch et al. 



errors were determined by generating Monte - Carlo data 1998). Finally, the A401 temperature profile that we use in 



sets which properly account for the statistical and system- 
atic uncertainties (including PSF effects). 

Figure ^ shows the best fit projected temperatures and 
la errors in each of the four radial bins. We note a slight 
(by a factor of 1.2 or 1 a) increase of the temperature in 
the center, compared to a quite constant value in bins 2'- 
5'-9'. The temperature in the 9'-16' bin falls by a factor of 
1.8 or 2.7 (T below the value in the 5'-9' annulus. Similar 
radial temperature decline s have been observed in a large 
sample of ASCA clusters ( [Markevitch et al. 199'^ ). The 
single temperat ure fit gives a mean tem perature of kT = 
8.0 ± 0.4 keV ([Markevitch et al. 1998|), which is consis- 



(1998 ) method (see Donnelly et al. 1999, Donnelly et al 



tent with the EINSTEIN MPC value kT = 7.8Io;g keV 



(David et al. 1993). The temperature profile values are 
consistent with the single temperature only by a proba- 
bility of smaller than 10"®, therefore A401 is significantly 
non-isothermal. 

We note that even though at the cluster center we saw 
a slight brightness excess (~ 30 %) compared to the (3- 
model, the sky area covered by the central bins (r < 0.5') 
is only 6% of the area covered by the innermost ASCA 
temperature bin of r < 2'. Therefore the contribution of 
the brightness excess to the total emission from the central 
temperature bin is negligible, and we ignore the effect of 
a possibly different temperature for this excess emission. 

The derivation of cluster temperature profiles from the 
ASCA data is not straightforward due to the wide and 
energy-dependent PSF. Takahashi et al. (1995) showed 
that if the PSF effect is neglected, an intrinsically isother- 
mal cluster will appear significantly hotter with ASCA at 
large radius (a 7 keV cluster would appear to have a tem- 
perature of ^ 20 keV at a 20' radius). For A401, Fujita et 
al. (1996) analyzed the same ASCA dataset and derived an 



this wor k, is similar to profiles of a large sample of ASCA 
clusters (Markevitch et al. 1998), when compared in the 
physically meaningful units of the radius of a given over- 
density. Therefore it appears unlikely that the observed 
temperature decline is caused by any unknown detector- 
dependent systematic effect. 

3.2. ROSAT 

In addition to our four ASCA temperature points, 
we have some crude temperature information from the 
ROSAT data. Since we have a significant (~ 3 a) detec- 
tion of the cluster gas in the 22'-25' (= 2.5-2.9 h^g Mpc) 
annulus we know that the gas temperature there must ex- 
ceed zero. Therefore we introduce a ROSAT temperature 
point kT > keV at 22' - 25'. 

4. MASS CALCULATION 

We will now use the ASCA and ROSAT temperature 
data to estimate the total mass of A401. For this, we 
assume that A401 is spherically symmetric and that its 
gas is in hydrostatic equilibrium (as indicated in Section 



2.1). From this condition the tot al mass within a sphere 



of radius r can be written as (e.g. Sarazin 1988) 



1-1 T(r) r ( dlnpnn, dhiT 

Mtot{< r) = 3.70xlO"Afo — ^-^ ^-^^ 

^- ' °10keVMpcV rflnr dlnr 

(4) 
using /I — 0.60. 

We consider the total mass consisting of stellar mass 
in galaxies, intracluster gas, and dark matter. We esti- 
mate the amount of stellar mass in galaxies using Dressler's 
(1978) King profile fit to the galaxy distribution in A401 
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with a core radius of 0.4 h^Q Mpc and a V band luminos- 
ity of 1.5 X 10^^ h^Q L0 inside 1 core radius. Even though 
the fit extends only to a 12' radius, we extrapolate this 
distribution to 15.1' {— rsoo) using the above King profile, 
and obtain a V band luminosity of 1.0 x lO^'^ h^Q Lq. We 
convert this value to galaxy mass, u sing a mass-to-ligh t 
ratio of 3.2 Mq/Lq /igo (taken from [White et al 19931 ), 
assuming a Coma-like luminosity function for galaxies in 
A401, and an M/L relation from van der Marel (1991). 
We find a stellar mass in galaxies of 3.3 x lO^'^ h^Q Mq, or 

16 /i^o % of our gas mass value, or 4 % of our total mass 
value within rsoo- Since this is much less than the uncer- 
tainty of our total mass estimates, we will not model this 
component separately, but rather include it in the dark 
matter model. This does not introduce any ambiguity in 
the interpretation of dark matter parameters (except for 
the normalization), since our best fit dark matter profile 
scales like the King profile at large radii (see below) and 
the core radius of this model is similar to the best fit King 
profile of the galaxy distribution in A401 (Dressier 1978). 
Quantitatively, the total mass is given as 



Mtoti< r) 



■i-n-r^ipdark + Pgas)dr, 



(5) 



where pdark and pgas denote the density profiles of the 
dark matter plus stars and gas, respectively. 

In principle, the mass can be calculated directly from 
Equation m if the gas temperature profile is known in de- 
tail. However, our temperature profile is not of sufficient 
quality to allow this procedure. We therefore use an in- 
direct method which assumes various models for the total 
mass radial distribution, calculates the corresponding tem- 
perature profiles, and compares them to the data, looking 
for acceptable models. Following Markevitch & Vikhlinin 
(1997b), we model the dark matter density distribution 
Pd using two different functional forms, which together 
approximate a wide range of physically motivated spheri- 
cally symmetric distributions. A constant core model has 
a dark matter density given by 



Pdark ^ 



1 



-a/2 



and the central cusp profile is described as: 



Pdark OC 



ad 



r 
ad 



rj — a 



(6) 



(7) 



In the equations above, ad is the scale length of the dark 
matter distribution. In the cusp model (Eq. |^), the first 
term describes the cusp behaviour near the center. In both 
models, the density at large radii scales as r~". With -q = 
1 and a = 3, the cusp model corresponds to the "universal 
density profile" which Navarro et al. (1995, 1997, NEW 
hereafter) show to be a good description of cluster CDM 
halos in simulations of hierarchical clustering. Since the 
quality of the data is not adequate for deriving the values 
of all the parameters independently, we fix ?/ = 1 in the 
cusp models, as suggested by the NEW simulations, but 
vary the other parameters. 

As in Markevitch & Vikhlinin (1997b), we solve the hy- 
drostatic equilibrium equation (Eq. H) for the gas tem- 
perature as a function of radius and gas and dark matter 



density parameters (using the (3 model for the gas density 
profile) : 



r(r)-(l+x2)W2 
where x = r/ax, 



To -A (l + y2)-3'3/2 



Hy) 



dy 



(8) 



A = 



All 



Pgasi^) 



3.70 X 1013 VMpcy Mq Mpc 



33keV (9) 



and 



Jo Pgas(0) Jn 



fd{z)dz, (10) 



where fd = Pdark{x)/pdi, and pdi is the dark matter den- 
sity at the X-ray core radius a^, and pdark{x) is given by 
either Equation (||) or (|^). Tq and Pgas(O) are the gas 
temperature and density at a; = . 

Eor the gas density distribution, we use the /3 model 
parameters deriv ed f rom the ROSAT surface brightness 
analysis (Section ^.2[ ). The remaining parameters to be 
fitted to the temperature data are a, ad, Pdi and Tq. 

While computing the temperature profiles, we use an- 
alytic solutions for the integrals in Equation ( pX)| ) in the 
cluster center to avoid the numeric effects of the singular- 
ity in the integral in Equation (||). Beyond 0.2 a^ we switch 
to numerical integration, preserving continuity. Step sizes 
for the numerical integrations are chosen to achieve a 1% 
accuracy in the computed temperatures. 

For each set of dark matter profile parameters and Tq we 
computed the 3D temperature profile, and projected it to 
the observed 2D ASCA annuli, weighting the line of sight 
temperatures with the emission measure of each volume 
element intersected by the ASCA annuli. The projected 
model temperatures were then compared with the mea- 
sured ASCA temperature values. The parameter values 
were changed iteratively to minimize x^. 

To incorporate the ROSAT temperature point (Section 
3.2) in the fitting procedure, we do the following: if the 
model temperature at 22' - 25' becomes negative, an ex- 
ponential increment is added to x^ (the more negative the 
value, the higher the increment), but as soon as the trial 
value in the fit becomes positive, the contribution of that 
data point to x^ vanishes. This arrangement will "steer" 
the fitting procedure smoothly towards positive tempera- 
ture values in the 22'-25' bin, without bias towards any 
arbitrary temperature value. 

The model temperature profile is very sensitive to small 
changes in the parameter values. On the other hand, the 
temperature errors and the widths of the radial bins are 
quite large, due to poor statistics. Therefore, a large range 
of parameter values gives a good description of the data. 
These two features led to difficulties using standard min- 
imization routines, which often find local x^-minima. We 
dealt with th is problem by app lying the simulated anneal- 
ing method (Press et al. 1994). 
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5. RESULTS 

Before presenting the results of our mass fits, we briefly 
address the question of convective stability, since some of 
our model temperature profiles, formally allowed by the 
data, have strong gradients at large radii and at the very 
center (especially the cusp models). These strong temper- 
ature gradients may not be consistent with the require- 
ment of convective stability. If the gas at radius r is to be 
stable against convection, the effective polytropic index at 
that radius, defined as 



dlogp(r) 



(11) 



must be less than |. During a dynamical time, all con- 
vective instabilities should have been erased, and clusters 
in general should be convectively stable, therefore strong 
temperature gradients cannot exist. However, at large 
radii, hydrostatic equi hbrium may not hav e time to es- 
tablish. Simulations ( [Evrard et al. 1996| ) suggest that 
''500 (the radius where the mean interior density is 500 
times the critical density) provides a conservative upper 
limit for the radius inside which the gas should be hydro- 
static. However, in many simulated clusters, hydrostatic 
balance holds at evenlarger radii. For A401, rsoo = 15.1', 
and ri5o — 25', the radius where the mean interior density 
is 150 times the critical density. Hence, at the radial range 
of our ASCA data (r < 16') the hydrostatic equilibrium 
assumption is likely to be valid, and we are justified to re- 
quire that the model temperature profiles be convectively 
stable. At the maximum radius of ROSAT detection (25') 
the case is less certain. In the Evrard et al. (1996) simu- 
lations, some clusters exhibit significant gas bulk motion 
at these radii. Depending on the cosmological scenario, at 
ri5o the kinetic pressure may be comparable to the ther- 
mal pressure, and the systems may be far from hydrostatic 
equilibrium. Since A401 is an apparently relaxed cluster, 
strong deviations from equilibrium are unlikely even at 
that radius. In order to be conservative, we will not apply 
the above convective stability constraint, 7 < |, at radii 
beyond r5oo. Also, we do not apply it to our solutions at 
r < 0.5' because at that radial range the cluster brightness 
does not agree with the (3 model. 

While using the ROSAT temperature information at 25' 
in our fits, we are extrapolating the hydrostatic tempera- 
ture profile to a possibly non- hydrostatic region of the clus- 
ter (out to ri5o), but this should not introduce a significant 
error. Cluster formation simulations show that tempera- 
ture profiles of clusters decline with the radius without 
any dramatic change in the temperature profiles between 
rsoo and ri5o, so that our extrapolation is justified, even 
though the kinetic pressure may become comparable to 
the thermal pressure at riso. Furthermore, if the gravity 
of the total mass inside ri5o balances the thermal plus ki- 
netic pressure at that radius, for a given mass distribution 
the temperature implied by the hydrostatic model is higher 
than the observed one. Hence, in our fitting the require- 
ment that the hydrostatic model temperature exceed zero 
at ri5o, is conservative. A possi ble complication, a Tp - T,; 
nonequality at large radii (e. g. |Fox & Loeb 1997|, lEttori & 



5.1. Core model 

We now discuss results of fits using the dark mass pro- 
file of the core model (Eq. g). To avoid local minima, 
we adopted a scheme to fix a at several values over an 
interesting range and to fit the other parameters. This 
was necessary, since the effect of this parameter is most 
significant at large radii where we do not have ASCA tem- 
perature data. The best fit core model is acceptable, with 
Y^ = 1.29 for 4 parameters and 5 data points (see Figures 
m and H and Table [|). The best fit core model gives a 
value for the dark mass M^ — 0.73 xlO"^^ H^q M© within 
r < rsoo, with Tq = 10.1 keV, pdi = 4.3 xlQ-^e g cm'^ 
(which corresponds to the central dark matter density of 
1.5 X 10-25 g cm-3), ad = 2.29 arcmin (= 260 h^f^ kpc), 
a = 3.1. The central dark matter density is 14 times that 
of the gas. The core radii of gas and dark matter mod- 
els are quite similar, and the dark matter density at large 
radii falls faster, scaling as r"'^-^, whereas the gas density 
falls as r~^-^. To check whether the solution agrees with 
the requirement of convective stability, we computed the 
effective polytropic index as a function of radius, using 
Equation O and used the gas distribution from the best 
fit P modeL This best-fit model is convectively stable in 
the r = 0.5' — 15.1' range. The overdensity in this best-fit 
model drops below 500 at radius 



rsoo - 15.1' = 1.7 h^^ Mpc 



(12) 



Fabian 199^ , Takizawa 199E) would have the same effect 



To propagate the errors of the temperature profile data 
to our mass values, we fit a large number of Monte - Carlo 
tem pera ture profiles with added random errors (see Sec- 
tion 3.1) using the same approach as for the best-fit tem- 
perature profile. We noticed that most models had a ten- 
dency for the temperature to fall to zero at radii below 25' 
(because of the strong decline of the temperature between 
the two outermost ASCA bins), therefore our ROSAT tem- 
perature information at large radius provides a powerful 
constraint. We rejected unphysical models which gave in- 
finite temperatures at large radii. Each Monte - Carlo 
model was checked against the convective stability con- 
straint, and rejected if 7 > | in the r = 0.5' — 15.1' range. 
The excluded models are those with the sharpest temper- 
ature peaks in the center, or those with the smallest core 
radii, which are also the models with smallest a due to the 
parameter correlations. The mass in the models with the 
smallest a values increases fastest with radius, therefore 
large masses are constrained more strongly. From the dis- 
tribution of the accepted Monte - Carlo models, we deter- 
mine the 1 a scatter of the mass values as a function of the 
radius. We convert these values to 90% confidence values, 
assuming a Gaussian probability distribution. Although 
in general we cannot constrain the dark matter model pa- 
rameters independently, some of these parameters are cor- 
related and the corresponding mass values vary within a 
relatively narrow range and can be reasonably constrained. 

Since the total mass is proportional to the local loga- 
rithmic derivative of the gas density distribution (oc (3 in 
Eq. (2)), we quantified the effects of the uncertainty of the 
local P value and its possible deviations from the global 
value (see Section |2.2D that we used in our mass computa- 
tion. We divided the ROSAT profile into radial ranges r = 
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0.5-5', 5-10', 10-20', 15-30' and fit these profiles with the 
(3 model, fixing the core radius to its global value. Within 
15' the P values are similar to the global value, while be- 
yond 15' the profile becomes slightly (but not significantly) 
steeper. We added the local /3 uncertainty in quadrature 
to the total mass uncertainty, which gave a very small ef- 
fect within rsoQ. Our final mass values for the core model 
at the gas core radius r = Ux (= 2.6' = 290 h'^^ kpc) and 
at r5oo (= 15.1' = 1.7 Mpc) are 



Mtoti< a.) = 1.20t°i^ X 10"/i5o' Mq 



(13) 



and 



Mtot{< rsoo) = 0.94t0;24 x lO^^/igV Mq. (14) 

Figure pi shows the resulting mass profile and the corre- 
sponding igas profile. 

5.2. Cusp model 

We now apply the cusp model given in Eq. 0. The best 
fit cusp model (77 = 1) has an unreasonably high value 
of the slope parameter a, due to parameter correlations. 
However, for the best fit models, the total mass within 
rsoo does not depend significantly on a. The best fit total 
mass within r^oo is 0.96 xlO^^ h'^Q Mq, almost identical 
to the best fit value for the core model obtained above 
(Eq. 1^. In the radial range a^ < r < r^oo the enclosed 
masses in the best fit core and cusp models differ by not 
more than 5%, a deviation that is negligible compared to 
the mass uncertainties. The scatter of mass values for the 
cusp model is smaller than that for the core model, so that 
at each radius the 1 a interval of the cusp model masses 
lies within that of the core model. Therefore, the mass 
values and errors obtained earlier with the core model will 
be our final values (see Table g). 

Even though the cusp model gives an acceptable fit to 
the temperature data, because of the centrally peaked form 
of this model it always violates the convective stability con- 
straint at the center, even at radii r > 0.5' where the gas 
profile is well defined by the (3 model. Models with 77 < 1 
have a less prominent peak, but as 77 approaches zero, the 
models essentially approach the constant core model al- 
ready considered in Section 5.1. The masses in the best 



fit models with different 77 values within rsoo equal that 
of the 7/ = 1 model. Since the gas distribution in A401 is 
well represented by a /? model inwards to a rather small 
radius (0.5'), one way to reconcile the cusp model, such 
as those predicted by simulations, with the observed gas 
density profile is to in troduce a significa n t non-thermal 
pressure in the center (Loeb & Mao 1994, Markevitch & 
Vikhlinin 1997). Such a pressure must satisfy the hydro- 



static equilibrium condition without inducing turbulence, 
which would require, for example, an equation of state 
Ptherm + Pnontherm oc p'^ with 7 < 5/3 (e.g., Landau & 
Lifshitz 1959). 

Note that A401 is rather unusual in that it has a cD 
galaxy but no significantly detected cooling flow and gas 
density peak usually found in cD clusters. However, the 
results of Peres et al. (1998) allow an upper hmit of 120 
MQ/yr for the mass flow rate in A401. Significant mass 
deposition from the cooling flow would make the average 
gas temperature less peaked, and possibly the NEW profile 
acceptable. Unfortunately, the quality of the current data 



is not adequate to construct a proper two-phase modelling 
of the cluster medium to test this possibility. 

5.3. Comparing m,ass values 

Eor comparison, we have calculated the mass profile 
under the traditional assumption of isothermality. The 
emission-weighted temperature model of the gas exclud- 
ing the co ntaminating component s for A401 gives Tjjf = 
8.0 ± 0.4 ( Markevitch c_^ al. 1998 ). Assuming a constant 
temperature. Equation W reduces to 



Mtot{< r) = 1.11 X 1015/3 



Tx 



{r/axf 



10 keV Mpc 1 + {r/ox 



-Mc 



(15) 

Figure pb shows this mass profile together with that de- 
rived using the observed temperature profile. At r = a^, 
(= 2.6' = 0.29 h^Q Mpc) the mass derived using the ob- 
served temperature profile exceeds the isothermal mass by 
a factor of 1.3, both agree at a radius of 11' = 1.2 h^^ Mpc 
and at 15.1' (= 1.7 Mpc = r^m) our value falls to 0.9 of 
the isothermal value. Qualitatively similar behaviour was 
found by Markevitch & Vikhlinin (1997) for A2256, where 
the effect had a larger magnitude due to the somewhat 
steeper temperature decline and wider radial coverage. 

Using the deprojection technique, and the broad beam 
EINSTEIN M FC temperature of kT = 7.8±0.6 keV ( |Da^ 
et al. 1993| ), White & Fabian (1995) and White et al. 
1997) obtained gravitating mass values of 1.01, 0.59 and 
1.07 X 10^5 /i^o^Mq inside radii of 1.265, 0.862 and 1.380 
hgQ^ Mpc, respectively. These masses are higher than 
the isothermal masses obtained by evaluating Equation 
(|l5|). Our temperature profile masses at the same radii 
are 0.74t^;J^, 0.53J:^;J!^ and 0.79l^;i? x lO^^ h^^ , respec- 
tively, which are lower by factors of 1.4, 1.1 and 1.3 (2.6, 
0.7 and 2.5 a). This difference is in line with our finding 
above that, with increasing radius, the isothermal model 
overestimates the mass. 

It is useful to compare our measurement with the cosmo- 
logical simulations of Evrard et al. (1996), who obtained 
cluster mass-temperature and radius-temperature scaling 
laws of the form: 

r5oo(7x) = 2.48 ±0.28 {Tx/lQkeVY/^ Mpc (16) 

and 

M5oo(Tx) = 2.22±0.55 (Tx/lOfceF)^/^ xlO^^ Mq. (17) 

For a cluster with Tx = 8.0 ± 0.4 keV, Evrard et al. 
(1996) predict rgoo = (2.22 ± 0.16) Mpc and Afsoo = 
(1.59 ± 0.25) x W^Mq. Using our best fit core model, 
the corresponding values (Eq. |l^ and |l^) are by fac- 
tors of 1.3 and 1.7 (3.1 a) lower than the values pre- 
dicted by the scaling law, respectivel y. In the case of 
A2256 ( [Markevitch fc Vikhlinin 1997| ), a similar differ- 
ence was found. The isothermal model (Equation 03) gives 
-^^500 = 1.1 X 10^5 /igQ^, a factor 1.4 lower than the scal- 
ing law value. The difference is due to two effects: the 
simulated clusters have steeper gas density profiles and 
shallower temperature profiles than those observed. 

Buote & Canizares (1996) studied the ellipticity gra- 
dients of the ROSAT data for A401 and derived the to- 
tal mass distributions in A401. Their shape for the total 
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mass density {ptot oc r ■*) is consistent with ours, but the 
normahzation is very different. Their total mass values 
within radii 0.8, 1.6, and 2.4 h^^ Mpc are 1.47-1.70, 2.26- 
4.00 and 2.56-5.73 xIQ^^K^Mq, while our values at these 



'50 



radii are 0.46-0.61, 0.71-1.10 and 0.73-1.51 xIO^'^/i^o^Mq, 
smaller by factors of 2-4, 2-6 and 2-8, respectively. This 
behaviour for A401 is similar to what Buote & Canizares 
(1996) find for A2199, for which they obtain 6 times larger 
total masses within 0.8 h'^^ Mpc, compared to a (isother- 
mal) /3 model estimate. 

In the optical, the virial theorem ( Girardi et al. 1998 ) 
gives R„jr 

hO.92 



4.6 h^Q^ Mpc = 40.0 arcmin) and M^.^, 



2.74j;gg2 xlO^^/igg-^Mo, while our values extrapolated to 
that radius are 1.76 l^'^s x 10^^ ft-^Q^ which is consistent 
within 90% confidence errors, but note that the extrapo- 
lated values are very uncertain. 

5.4. Baryonic fraction, flm 

The best fit models show that the dark matter den- 
sity declines more rapidly than the gas density, which also 
means that the gas mass fraction fgas{< r) = Mgas{< 
r)/Mtot{< f) is a monotonically increasing function of ra- 
dius (see Figure ||c). In A401, at rsoo the gas mass fraction 
reaches a value of 



fgas{< rsoo) - 0.21 + 



+0.06 1-3/2 
-0.05 "SO ■ 



(18) 



This behaviour and value are similar to the results of a 



sample of EINSTEIN clust ers ([White fc Fabian 1995] ), a 
sample of ROSAT clusters ([DavMet al. 1995| ) and A2256 
( [Markevitch fc Vikhlinin 1997| ). As shown by White et 
al. (1993), fgas has important implication for the cosmo- 
logical matter density parameter f2„i —< p > / Pcrit- We 
define T as the ratio of the local baryon fraction /h in a 
cluster to the primordial value r2h/f2„i. Therefore, 






-T = h- 



(19) 



Assuming that the baryonic matter consists of the gas and 
stellar mass in the cluster, we write 



h 



Mtnt 



Hence, 



^ 't71 J^ ^ '6 [ J gas 



J gas + 



Mgai 
Mtnt 



Mgai 

Mtnt ■ 



(20) 



(21) 



We evaluate Equation ^ at rsoo using 1) our fgas value 
(Equation p^, 2) our galaxy mass estimate from Section |[ 
3) a slight baryon diminution T(500) — 0.90 suggested by 
simulations, (Frcnk et al. 1996), 4) f}b/i|o = 0.076 ±0.007 



(Buries et al. 1998), and propagate the errors of fgas and 



VLb- Figure [fl shows the resulting allowed parameter space 
of (fimjilo)- If some dark matter is baryonic, then J7,„ 
would decrease further. fi™ = 1 is allowed only by an un- 
realistically low value of the Hubble constant, Hq < 8 
km s~^ Mpc^^. Using, for example a value Hq = 68 ± 8 
km s~^ Mpc^^ found from the analysis by Nevalainen & 
Roos (1998) who studied the Cepheid metallicity effect on 



galaxy PL- relation distances calibrated at LMC, we find a 
cosmologigal matter density parameter (90 % confidence) 
of 



which is consistent with the value obtained by combin- 
ing all relev ant current independent ^rn est imates, ^m = 



0.33 ±0.11 ([Roos fc Harun-or-Rashid 1999[) 



6. CONCLUSIONS 

Using spatially resolved ASCA spectroscopic data, we 
have constrained the dark matter distribution in A401, 
without the assumption of isothermality. The dark mat- 
ter density in the best fit "constant core" model scales 
as r~^-^ at large radii. Thus a well-known King model 
appears to describe the dark matter distribution, as well 
as the galaxy distribution in A401 (Dressier 1978). This 
slope is also the same as found by NWF in their simula- 
tions, although the simulated clusters also exhibit a cen- 
tral density cusp. For A401, such a cusp violates the con- 
vective stability condition in the cluster center, because 
the gas density is well described by a standard /3 model. 
One way to reconcile a total mass profile having a cusp 
shape with the observed gas density profile is to introduce 
a significant non-thermal pressure in the center. Such a 
pressure must satisfy the hydrostatic equilibrium condi- 
tion without inducing turbulence, which would require, 
for example, an equation of state pthcrm + Pnonthcrm oc p^' 
with 7 <C 5/3 (e.g.. Landau & Lifshitz 1959). Alternately, 
significant mass deposition from the cooling fiow would 
make the temperature less peaked, and NEW profile ac- 
ceptable. However, the quality of the data is not adequate 
to test this hypothesis. Regardless of the presence or ab- 
sence of a central cusp, the total mass within rsoo (1-7 /i^o 
Mpc) is 0.94l°;^4 x lO^^ h^^ Mq at the 90 % confidence. 
The mass within the X-ray core (290 h^^^ kpc) is a factor 
of 1.3 higher than the value from an isothermal analysis 
while at rsoo the value is 0.9 the isothermal value, which 
is qualitatively similar to the A2256 result (Markevitch 
& Vikhlinin 1997). Our Msoo value is 1.7 times smaller 
than that predicted by the scaling law of Evrard et al. 
(1996). This discrepancy arises because the simulations 
do not correctly predict the observed gas density and tem- 
perature profiles. The gas density profile is shallower than 
that of the dark matter, being proportional to r~^-^ at 
large radii. Hence the gas mass fraction increases with 

radius, with fgas{r-,m) = 0.21+o:o5 ^ot ^^ (90 % errors) at 
''soo- Assuming that the cluster matter content is repre- 
sentative of that of the Universe, this implies Q,m < 0.31 
at 90% confidence, in conflict with the Einstein-deSitter 
Universe. 
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Table 1 
Best-fit values with 90 % confidence errors 



Parameter 


Value 


Qx [arcmin] 


2.56 ±0.14 


Ox IKfl kpc] 

/3 


294 ±16 
0.70 ±0.02 


Iq [10^^ c s^-^ arcmin ^] 
xVd.o.f. 


5.6 ±0.3 
49.7 / 55 


PgasiO) [1014 hli'' Mrs MpC-3] 
P,as{Q) [lO-^*^ hli' g Cm-3] 


1.6 
1.1 


To [keV] 

pM [10-25 g cm-3] 

ad [arcmin] 


10.1 

1.5 

2.29 


ad [Ko kpc] 
a 


260 
3.1 



'using Hq = 50 /i5o km s ^ Mpc 
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Table 2 
Mass values with 90 % confidence errors 



Radius Mgal Mgas Mtot fgas X h%^ 

[/igo'Mpc] [1013 /jg-i Mo] [IOI^/j-IMq] [1015/1-1 Mo] 



0.29 (= a,) 
1.0 

1.7 (= 7-500 ) 




0.23 
2.0 
3.3 


0.11 
0.99 
2.01 


n 120+°°" 

U.iZU_o.o49 

oMtrot 


0.09l°o;St 

o.i6t°o;^^ 
o.2itS;;!^ 


* using Hn = 


;50 /i5 


km s^i ]\ 


/Ipc-i 
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Fig. 1. — Contour map of the surface brightness of combined rp800182n00 + rp800235n00 pointings, smoothed by a Gaussian with a = 1'. 
The contour level values are 0.0005, 0.001, 0.002, 0.004, 0.008, 0.016 and 0.032 counts s^^ arcmin"^ ^401 is located in the center of the 
image, and the neighboring cluster A399 is located south-west from the image center. 

Fig. 2. — ROSAT PSPC radial surface brightness profile together with the PSF-convolvcd best-fit /3 model. The data with r < 0.5' are not 
included in the fit. Note the slight excess in the center. 

Fig. 3. — Confidence contours for the a^ and (3 parameters of the surface brightness profile fit in Figure 2. 

Fig. 4. — Crosses show projected ASCA temperatures with 1 cr errors, and the lower limit from the ROSAT data. Thin solid lines show a 
representative set of temperature models (before projection), allowed by the convective stability constraint (see text). Thick solid line shows 
the best fit model (before projection). Thick dotted line shows the values of this best model projected to the 2D ASCA bins, which are 
compared with the ASCA data. 

Fig. 5. — Mass distributions of the models plotted in Figure 4: (a) the dark matter and the gas densities, (b) the enclosed total masses 
with the errors and the gas mass and (c) the gas mass fraction with errors. For comparison, values assuming isothermality (kT = 8.0 keV) 
are also shown. Masses are evaluated using Ho = 50 km"-*^ Mpc~^ (total mass scales as H~^ and gas mass as H~'''^). 

Fig. 6. — The 90 % confidence area in parameter space of [Ho,Qm], using the derived value of fgas- A very low value of Ho is needed for 
Qm= 1. The shaded area shows the subspace allowed by Ho value from Nevalainen & Roos (1998) 
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A401 data and folded model 
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